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Abstract 

We propose a method to construct the ground state ip{X) of local lattice hamil- 
tonians with the generic form Ho + XHi, where A is a coupling constant and Hq is 
a hamiltonian with a non degenerate ground state V'o- The method is based on the 
choice of an exponential ansatz ^(A) = exp(U{X))'ipo, which is a sort of generalized 
lattice version of a Jastrow wave function. We combine perturbative and variational 
techniques to get succesive approximations of the operator U{X). Perturbation the- 
ory is used to set up a variational method which in turn produces non perturbative 
results. The computation with this kind of ansatzs leads to associate to the original 
quantum mechanical problem a statistical mechanical system defined in the same 
spatial dimension. In some cases these statistical mechanical systems turn out to be 
integrable, which allow us to obtain exact upper bounds to the energy. The general 
ideas of our method are illustrated in the example of the Ising model in a transverse 
field. 







1 Introduction 



In Statistical Mechanics and Quantum Field Theory there is a large class of two-dimen- 
sional models whose basic magnitudes can be computed exactly|l[|. The mathematics and 
physics of these integrable models are very rich, constituing areas of intensive investiga- 
tion. There are in contrast very few integrable models in higher dimensions^, ^, where 
almost all the models which are interesting from a physical point of view are believed to 
be non integrable. These facts represent a major conceptual and technical gap between 
the 2d and the higher dimensional worlds, which would be desirable to fill in. 

The lack of exact techniques has always motivated the construction of approximative 
methods which very often give good quantitative and qualitative results. These methods 
can be roughly classified according to four main categories, namely: perturbative, varia- 
tional, numerical and renormalization group, every one having advantages and limitations. 
The perturbative methods can be developed in a systematic way but are restricted to small 
values of the coupling constants. The variational ones give in general non perturbative 
results but rely on aprioristic or intuitive conjectures. The numerical approaches are in 
general unbiased but limited to small lattices. Finally, the RG techniques when combined 
with one of the previous ones are very powerful for the study of critical phenomena. 

In this paper we shall restrict ourselves to the first two methods. Our aim is to combine 
perturbative and variational techniques choosing the best of each and by passing some of 
their limitations. The problem we cope with is the construction of the ground state of 
a local lattice hamiltonian of the form Hq + XHi, where A is a coupling constant. Our 
strategy is to use perturbation theory to set up a variational approach to the computation 
of the ground state wave function. The starting point is the choice of an exponential ansatz 
for the exact ground state of the hamiltonian Hq + XHi. The peculiarity of this ansatz is 
that, under certain assumptions about the ground state of the unperturbed hamiltonian 
Hq , the whole hierarchy of perturbative equations can be solved systematically in terms 
of a collection of irreducible local operators Vj, which are characterized by the order uj at 
which the operator Vj first appears in perturbation theory. Each operator Vj comes in the 
ansatz multiplied by a weight q;/(A) = higher powers. The data {Vj, a/(A)} contain 
all the information required to built up the ground state. For small lattices it is easy to 
construct this set and therefore obtain the exact ground state, but for large lattices we 
have to resort to some kind of approximation. The natural thing to do is to truncate the 
ansatz considering only those operators Vj with a level equal or less that a given number 
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u, which is the order in perturbation theory where one is working at. The problem is 
then to find the weights a/ (A) for z// < u. At this point one can follow two different 
paths: i) use the truncated ansatz as a trial state in a standard variational fashion or ii) 
take the truncated ansatz as the exact ground state of a hamiltonian which should be 
sufficiently close to the original one. In this paper we shall develop the approach i) and 
leave the approach ii) for another publication. We have so far outlined the main ideas 
and techniques which will be exemplified in the case of the Ising model in a transverse 
magnetic field (ITF), but it should be clear from the exposition that they apply to more 
complicated models as for example the chiral Potts model, the XXZ- model, etc. 

Trial wave functions with an exponential or, more generally, product structure are 
known as Jastrow wave functions and has been applied to a wide variety of problems in 
different areas as for example, liquid ''He 0] , nuclear matter P , solid state physics p, 0] 
and more recently in the fractional quantum Hall effect [§. In quantum mechanics a Jas- 
trow wave function is simply the exponential of the potential energy, ip{x) = exp{—aV{x)) 
and becomes exact for the simple harmonic oscillator. In a many body problem one re- 
places V{x) in the exponential by J^i.jfi^i ~ ^j) where / is some two body effective 
potential. In this paper we shall consider discrete versions of generalized Jastrow states 
in contrast of the continuous versions described above. To our knowledge there is not in 
the literature a general and systematic study of this type of wave functions except for 
some interesting observations and discussions confined to particular models 0. Among 
these observations is the fact that the norm of the Jastrow wave functions coincides, in 
certain cases, with the partition function of an associated classical statistical mechanical 
model defined in the same dimension, with the variational parameter playing the role 
of inverse temperature. Indeed, in the example of the ITF we find that the norm of 
the exponential ansatz is given by the partition function of a classical Ising model with 
couplings among the spins dictated by the familly of operators Vj and as many coupling 
constants as weights aj. Working with these generalized Jastrow wave functions requires 
the knowledge of thermodynamic quantities such as internal energies, susceptibilities, etc. 
The calculation of these quantities is in general a difficult task and this imposes limita- 
tions to the variational approach i) mentioned above. However in certain occasions the 
computation can be carried out exactly. This happens when the associated statistical 
system which underlies a particular ansatz is integrable. If the quantum hamiltonian is 
one dimensional then the statistical system is also one dimensional and the corresponding 
partition function can be easily computed. Things are more interesting in two dimensions 
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where, as we said in the introduction, most of the hamiltonians are not exactly solvable. 
The approximation method that we propose bring us to statistical systems in the same 
dimension, and therefore in two dimensions we may some times make contact with 2d 
integrable models. This opens the possibility of using the huge amount of information 
known in integrable 2d systems in the study of 2d quantum systems. 

The paper is organized as follows. In section 2 we introduce the P/V method. In 
section 3 we apply the general ideas presented in section 2 to the Ising model in a transverse 
magnetic field, making special emphasis in the Id and 2d cases. In section 4 we present 
our conclusions and future prospects. All the technical details and computations are 
included in four appendices. 

2 The Pert urbative/ Variational Method 

Let H{X) = Hq + XHi be a local lattice hamiltonian where Hq is a Hamiltonian whose 
ground state ipo is non degenerate. We suppose that ipo and its energy Eq are known 
exactly. If A is sufficiently small one can use standard perturbation theory to expand the 
ground state wave function ipi^X) and energy E{X) of H{X) = Hq + XHi as follows. 




where {ipn} and are subjected to satisfy a set of equations linking quantities of the 

same order ||T0[. In order to write this set of equations we shall assume that ip{X) can be 
reached from ipo with the action of an operator of exponential type. 



^(A) = e^W ^0 



n 



(3) 



The relation between the perturbative expansion (|I|) and the operators Un is given by 
a non-abelian version of the Schur polynomials. 



^1 




(4) 
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V'n ~ Sp=l pi SjiiH \-np=n ^ni ' ' ' f^rip V'O 

If the operators f/„ commute among themselves then the above expressions coincide 
with the Schur polynomials Si, S2, Sn a^s functions of the variables f/„. The equations 
satisfied by ipn and E^"-^ translate now into a set of equations for Un and E^"'^ which we 
pass to derive. 

Introducing (0) into the eigenvalue equation, 

H{X) ^(A) = E{X) ^{X) (5) 
and multiplying on the left hand side by the inverse of the exponential one gets, 

exp £ A"f/„^ H{X) exp (^f] A"[/„j 7/>o = ^(A) t/^q (6) 
To simplify (^) we use the identity pT| , 



00 1 

e-^i?e^ = e^-(i?) = E-T^A(i?) (7) 



n=0 



where 



which yields, 



Ca{B) = [B,A], C\=[[B,A],A],... 



exp (^^ A"£c;„ j HiX) ^0 = ^^(A) V^o (9) 

The equations for the operators f/„ are finally obtained expanding the exponential in 
and collecting all the terms of the same power in A , 



Ho Vo = ^0 ^0 
{[Ho, U,] + i[[i/o, Ui], Ui] + [Hi, Ui]) ^0 = E('^ 



(10) 



J2p=lJ2ni-\ \-np=n pi ^U„-^ ' ' ' ^U„p{Ho) 

+ E;=iE„,+...+n,=n-l ^^C/„, ■■■Cu^^m] i,o = E^-)^o 



for n > 2 
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The solution of these equations is by no means unique. In particular the hermiticity 
properties of U{X) are not fixed a priori. There are two main choices that can be made in 
terms of hermitean or antihermitean operators. If U{X) is choosen to be antihermitean, 
then eq.(|]) becomes a unitary transformation which preserves obviously the norm of the 
state, i.e. < ipWl^piX) >=< V'olV'o >= 1- Moreover if that Hi has an off diagonal form, 
then E^^^ vanishes. When this happens it is customary to perform a unitary transforma- 
tion. 

Ho + XHi ^ e-'^^'iHo + XHi)e'^^ = H,s (11) 
where S' = S'"'" is obliged to satisfy, 

i[Ho,S]+Hi = (12) 

in order to cancel the terms proportial to A in H^s- This equation coincides with the first 
order equation in (|1^) but expressed in a matrix form, i.e. the dependence on the state 
ipo has been dropped. Finally the "effective hamiltonian" Hea has the expansion, 

H,f, = Ho + '-^[Hi,S]+OiX') (13) 

These kind of unitary transformations are extensively used in Statistical Mechanics and 
Condensed Matter Physics[jl2|, and have also been considered from a more mathematical 



point of view |jT3| 



In this paper we are interested in the study of those solutions to equations (|TD|) which 
are hermitean. An immediate consequence of this choice is that the norm of the state 
ip{X) is no longer the unity and in fact it will be an important quantity to compute, 

Z{X) =< >=< ^o|e^^(')|V^o > (14) 



A peculiarity of equations (|10|) is that all the operators [/„ appear always involved 
in nested commutators. Let us suppose that Hq and Hi are local operators defined on 
a lattice, i.e. operators involving a finite number of lattice variables located within a 
neighbourhood, whose size is independent of the size of the whole lattice. This kind 
of operators go, in the thermodynamic limit, to local operators in the continuum. We 
shall also assume that the number of degrees of freedom per site is finite. We want 
to investigate under which conditions the solution of the equations (|10[) can be given 
by a family of local operators. That this kind of solutions may exist is suggested by 
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the property mentioned above about the nested structure of equations (|ToD- The reason 
is that the commutator of local operators is again a local operator and therefore each 
equation (^) reduces to a sum of local operators acting on the unperturbed ground state. 
The other ingredient we shall require to achieve a local solution is an unperturbed ground 
state ipo with an ordered structure, which can be ferromagnetic, antiferromagnetic, or of 
some more general type. This is sufficient to guarantee that the operators Hq, Hi and 
Un should produce local disturbances acting on ipQ. There are many physical systems 
where these conditions are fullfiUed. The reason for requiring an ordered ground state 
is that in this case equations (p!0[) become a system of linear equations with a number 
of unknowns depending on the dimension of the lattice but not on its size. It will also 
depend on the order of perturbation theory and the number of lattice variables at each 
point. The aim of the above considerations is to make plausible the existence of solutions 
of the perturbative equations in terms of local operators f/„, but of course it does not 
constitute a proof. Later on we shall go through particular examples to show that this 
kind of solutions do indeed exist. In the case of the Ising model , and we suspect it is a 
general feature , every operator Un is given by a sum of irreducible local operators which 
we denote by V/, 

f/n = E Pl,n Vl (15) 

where the coefficients p/,„ are obtained by solving equations (p^O]). and where each operator 
Vj is a product of Pauli matrices located at different points of the lattice. The sum over 
/ in (|T5|) is finite and contains in general operators that have already appeared in the 
solution of equations of lower order. It is convenient to define the level uj of Vj as the 
order in perturbation theory at which this operator appears for the first time in the 
solution of equations (|lO|), 



z// = minimum {n / „ 7^ 0} (16) 

Introducing ([T5|) into (H) and interchanging the order of the sums J2n J2i = J2i J2ny 
we rewrite the exponential ansatz as follows, 

^{X) = exp (^^ E «/(A)^/) i^o (17) 
where ai{X) is the "weight" of Vj in the ground state. 
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aj{X) = 2 £ p,,„A" 



n=i/j 



The factor 1/2 in eq.(0) has been introduced for later convenience. Definition ([18 



imphes that a/ (A) has a Taylor expansion starting at X'^' , and containing non vanishing 
terms of higher powers whenever Vj appears in the solution of eqs.([TO|). The knowledge 
of a single function ai{X) will in general require solving the entire hierarchy ([l0|)!. At this 
stage one may wonder whether there are other alternatives to compute the functions «/ 
rather than solving the whole hierarchy (^^). To handle this problem it is illustrative to 
compare equations (|^) and ([T7|). They both give two different perturbative approximations 
to the ground state. While (^ is an expansion in the coupling constant A, eq. (p!7D is 
a kind of cluster expansion. In a sense they are conjugated perturbative expansions. 
Exploiting this analogy we shall propose an alternative way for computing a/ (A). In 
standard perturbation theory one truncates the state excluding corrections beyond a given 
power A". The analog of this for (|l^) will be to truncate the sum over I to contain only 
those operators Vj with a level uj less or equal than a given order i/, 

^('^)(A)=expfi Y: "i'^^W^/l A (19) 

In eq.(0) we leave open the possibility that the weights a^^'' depend on the order of 
the approximation u. The problem now is how to determine a^^\ In this paper we shall 
use a variational method which consist in taking (|19|) as a trial wave function for the 



ground state of the hamiltonian Hq + XHi. af^ are of course the variational parameters. 
Minimization of the energy of this state will give the weights af'^ as functions of A. 
It is clear that this minimization procedure should agree with the result obtained in 
perturbation theory to the same degree of accuracy, i.e. 

af\\) - ai{X) = 0(A"+^), V/ / z/, < z/ (20) 

On the other hand, it is well known that if one takes the perturbative solution to an 
order z/ as a trial wave function then the variational energy agrees with the perturbative 



expansion to order 2z/ + 1 |I4 



^(.ar,.)(^) _ _ 0(A2'^+2) (21) 
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This fact already illustrates at the perturbative level the usefulness of combining both 
perturbative and variational technique^. Next we consider the non-perturbative aspects 
of (|19[). The question is whether this ansatz can be extended to large values of A remaining 
close to the exact ground state. In particular in the limit A — > cxd we wish the state tp^'^'* 
to flow to the ground state manifold of the Hamiltonian Hi, i.e. 

lim Hiij^''\\) = El lim V^^^^H^) (22) 

A— >oo A— >oo 

where Ei is the ground state energy of Hi. If this equation holds then ip^'^\\) will 
interpolate in a continuous way the ground state of Hq and those of Hi, giving certain 
validity to the results obtained for the intermediate values of the coupling constant A, 
where a weak or a strong coupling expansion are dubious. A sufficient set of conditions 
for eq.(P^) to hold is the existence, within the family of operators expanded by {Vj}, of 
an hermitean operator J2i,ui<uo having a common eigenvector, say 0o, with the 

Hamiltonian Hi, and such that 0o and ipo have non- vanishing overlap, 

^) J2i,ui<iyo ^^1^°^ ^ is hermitean 

u) 00 is a ground state of Hi and Y,i,ui<uo (23) 
in) < V'ol^o >7^ 

Under these conditions it is easy to see that the ansatz for u >vq flows to 0o in 

the limit \ ^ oo. The interest of conditions (|23| ) is that they strongly restrict the possible 
solutions of the perturbative equations (|1^) (see appendix C for concrete examples). 



3 The Ising Model in a transverse field 
3.1 General Considerations 

In this section we shall study in detail the Ising model in a transverse field. Let us 
introduce our notations and some generalities valid in any dimension. 

The model is defined in a hypercubic lattice in d dimensions and L sites with a 



Hamiltonian given by |[T5| , 16, |T]^ 

^Recall that in perturbation theory the knowledge of the n^^ order correction of the wave function 
allows the computation of the energy to order n+1 according to the formula =< t/iqI^^iI^'-"'' > ■ 
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H - H f^fn 



(24) 



where and are Pauli matrices acting at the site n, with n = (rii, . . . ,nd 
1, . . . , L^^'^ and yU is any of the lattice vectors yUi = {1,0, ... ,0), fid = (0, 0, . . . , 1). 
The cases d = 1 and 2 correspond to a hnear chain and a rectangular planar lattice 
respectively. We shall suppose periodic boundary conditions. 

The phase diagram of this Hamiltonian as a function of A and the temperature is well 
known ||18| . We shall restrict ourselves to zero temperature, whose analysis is interesting 
because it shows a simple but non trivial example of quantum critical phenomena |T9|, 20| . 
In this case there is a critical transverse field Ac below which there are two ferromagnetic 
ground states characterized by the non vanishing of < >. For A > Ac the ground 
state is disordered and < >= but < >^ 0. Near Ac the critical behaviour 
of the \d Ising model with transverse field as a function of A at T = is the same as 
the critical behaviour of the classical Ising model in two dimensions as a function of T 
| T9| . This can be seen from the fact that the Ising model hamiltonian in Id is related to 
the transfer matrix of the 2d Statistical Ising model. The analogy between the quantum 
critical behaviour of the d-dimensional ITF and the statistical critical behaviour of the 
Ising model in d+1 dimensions for d=l is believed to hold for d > 1 ^ 22 1. This 



fact is not peculiarity of the ITF but a property which holds more generally |23, 23 1. 



In \d the model has been solved exactly by many different methods p6|, |T9l. The 
ground state energy of a periodic chain of L sites is , 



(Id) 



L/2-1 

E 

n=-L/2 



1 + A2 + 2A cos 



(2n + l)7r 

z 



which gives the following energy density in the thermodynamic limit, 



(25) 



^0 



lim 



(Id) 



' \ 2'2' ' (1 + A)2^ 



Eqs.(^5D and (^) exhibit the Krammers-Wannier duality. 



(26) 



eo ^(A)=Ae(^'^)(-) 



1, 



(27) 



in agreement with the fact that Ac = 1 is the selfdual critical transverse field. The density 
energy at the critical field is. 
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= -- (28) 
vr 

For dimensions higher than one the exact solution is not known and only approxima- 
tions are available In the next subsections we shall obtain variational upper bounds 



for the exact energy. It is also possible to obtain a lower bound for the ground state 
energy of the Hamiltonian (plD, 

4^^(A)>c^er(^) (29) 

There are two extreme cases where one can find the exact ground state of the Hamil- 
tonian (|2^. At A = there are two ferromagnetic ground states corresponding to 
< >= ±1 given by, 

I T>= ®n I T>n , I i>= ®n | i>n (30) 

where | !>„ , | i>n denote the eigenstates of with eigenvalues 1, —1 respectively. The 
density energy of these states is, 

eo(A = 0) = -d (31) 
At A = +00 the ground state is the non degenerate disordered state ( < cr^ >= ), 

|0 >= ®„|0 >„= ®n -^(1 t>n +1 i>n) (32) 

whose energy is given in the limit A ^ oo by, 

eo(A +00) ~ - A (33) 



The hamiltonian (24) commutes with the spin rotation operator, 

R = l[ (34) 

n 

and therefore its spectrum can be split into an even {R = 1) and odd {R = —1) sectors. 
The ground state ( P^D of the disorder phase belongs to the even sector whereas in the 
ordered phase there is a ground state belonging to each of the two sectors. 

We shall study the ordered and disordered phases independently. In the next two 
subsections we present the results of the computation of the energy and the variational 
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parameters using the variational method. The details of their derivation will be presented 
in the appendices. The formulae for the energy and the weights aj are labelled by the 
type of approximation used, its order and the spatial dimension of the lattice, namely, 
e(app, order, d), a/(app,order,d), where app= per (perturbative), var (variational). 

3.2 Ordered Phase 

In this phase the ground state is degenerate and therefore we could not in principle apply 
the techniques of section 2. However in practice this is not an obstacle since for large 
lattices the state | |> does not mix with the state | |> at finite orders in perturbation 
theory and one can effectively apply non degenerate perturbation theory, say to the state 
I t > [0 • In the ordered phase A is a small parameter and therefore the hamiltonian (p^) 
split as follows, 

H= Ho + \H^ 

Ho= -EnE^ (35) 

Hl= - En O'n 

The unperturbed ground state is taken to be, 

V^o = I T>=ni T>n (36) 

n 

A solution of the perturbative equations (0) which satisfy also conditions (p3D is given 

by, 

Ui= i.En^^ (37) 



^2- 16d2(2d-l) ^n,M^n ^n+M (2^) 

In appendix C we consider another solutions to equations ([ToD , which do not satisfy 
(p3|), and discuss their significance. 

The perturbative energy to third order in A which can be derived from this solution 




e(per,n = 3,d) = - + - (39) 
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Notice that odd powers in A do not enter into the energy. 

We shall now apply the P/V method corresponding to = 1 and 2. From eq.(p7D the 
ansantz for = 1 is given by, 



exp(|E„cr^)|T> 



(40) 



(cOShf I t>n +Sinh| | i>, 

which shows that ip^^^ is a mean field state parametrized in terms of hyperbolic functions 
of the variational parameter h. This mean field state is usually written in terms of an 
angle 6, which describes the deviation of the up pointing spin under the action of the 



transverse field |15|. These two kinds of parametrizations can be explained from the 
hermiticity properties of Ui, i.e. /;,(non-compact variable) ^ f/i(hermitean), ^(compact 
variable) ^ [/i(antihermitean) (see appendix C). 

The minimization of the energy of (^OD is achieved for, 

i -(d+^) X<2d 
e(var..= l.d) = j ^-^^ (41) 

For X < 2d the variational energy (|41|) coincides with the perturbative result (^9]), in 
agreement with condition (0). For X > 2d this energy coincides with the strong coupling 
limit (^), in agreement now with eg. (P^ . 

The mean field state (|40|) illustrates in the simplest possible way the interpolating 
properties of the variational states constructed out from solutions of the perturbative 
equations (|T0|) fulfilling conditions (|2^). Later on we shall see more examples of this. 

The dependence of the variational parameter on A is given by. 



arctanh (^^^ X < 2d 



h{y^T,u=l,d) = { V2dy - (42) 

oo X>2d 

The value = oo for A > 2(i implies that ip^^^ can be identified with the state |0 >, 
which is precisely the eigenstate of the Hamiltonian Hi. 
The magnetization of (HO) is. 



21 1/2 

A < 2d (^3j 

X>2d 



which shows that Ac = 2d and the exponent /?, defined by < >~ (Ac — A)^, takes the 
mean field value (3 = 1/2. The exact result in Id is given by. 
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^ J A<1 

<a^ >= { ^ ' - (44) 

[ A > 1 ^ ^ 

To go beyond the mean field approximation we shall consider the second order ansatz 



2, which according to (|38|) reads 



= (f E + f E ^n^-.. j I T> (45) 

The term in the exponential of (^) proportional to a introduces in t/''^^) correlations 
which were absent in '^'^^^ . Hence ■j/'''^-* is no longer a mean field state. The norm of ■j/''-^^ 
coincides with the partition function of a d-dimensional classical Ising model with inverse 
temperature a and magnetic field h. In \d the computations can de done exactly (see 
appendix A). The results for the energy and magnetization are shown in figure 1 and 2. 
The variational estimation of the energy for v = 1 considerably improves the mean field 
value result. The magnetization is also improved, but the singular mean field behaviour 
is lost obtaining instead a smooth curve (see table 1). This implies that in order to obtain 
an estimation of the exponent /? we should resort to another techniques as for example 
Fade approximants. The study of these questions, though interesting, is beyond the scope 
of this paper. 

In (i > 1 the norm of the state (^) cannot be computed exactly hence we have to use 
some kind of approximation. For small values of A we expect both h and a to be also 
small, in which case this norm can be computed through a "high temperature" expansion. 
Doing this computation we obtain, 

e(var.. = M) = -(d+^ + 5j^5^ + 0(A»)) (46) 

The fourth order term in A agrees with the perturbative result in accordance with (^Tj). 

For large values of A we expect that /i, a or both become also very large in which 
case the norm of the variational state can be computed by means of a " low temperature" 
expansion. Of course if the model is integrable we can compute in the whole range of 
couplings including the intermediate ones , obtaining rigorous upper bounds for the exact 
energy. 
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3.3 Disordered Phase 



In this phase the ground state is non degenerate and hence the techniques explained in 
section 2 could be automatically applied . The perturbative parameter in the disordered 
phase is j, which leads us to write the Hamiltonian (p4D as, 

Ho= -En (47) 

The unperturbed ground state is, 

V>o = |0>=niO>n (48) 

n 

A solution of the equations (p!OD and (|^) reads, 

Ui= iEn,, (49) 



The perturbative energy to third order is, 

e(per,n = 3,d) = -(A + -^) (51) 

4A 

In this phase we shall only study the first order ansatz (z/ = 1) which according to 
(1^) becomes, 

^(1)= exp(fEn,;.crfcrf_,J|0> (52) 

The norm of ijj^^^ coincides with the partition function of the d-dimensional classical Ising 
model with inverse temperature a. This implies that computations with (|5^) can only be 
done exactly in = 1 and 2. 

The Id case is very easy to handle ( see appendix B). The energy of (|52|) is after 
minimization, 

e(var, J. = l,ci =!) = <( " + ^) ^ " (53) 
^ ^ ^ -1 A<l/2 ^ ' 
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Compairing this result with the corresponding one in the ordered phase (|4l|) , we 
observe that they satisfy the duahty relation p7|), which means that we have incorporated 
the KW duality in our ansatzs. We have now two kinds of variational results at order 
z/ = 1, which can be combined in order to optimize the upper bound to the energy, 



l,d=l)= (54) 
Minimum (e(°'^'^'='^)(var, u = l,d= 1), e(<^''^°'^^^'^)(var, u = l,d = 




(55) 

A > 1 ^ ^ 

Let us define Aq]? as the crossing point between the ordered and disordered variational 
energies at order v. We expect that Aqd will approach Ac as z/ ^ cxd. Eq.(|^) shows that 
y^Qjy = 1 in Id. In fact from KW duality it is clear that Aq]) = 1 for all v in Id. 

Let us turn to the more interesting case of 2d, where the partition function and other 
related quantities are known exactly. In figure 3 we show our results for the energy. The 
value of Aq£) is 2.7239 which is not far from the critical value computed by other methods 



which is around 3.08 20 . 



3.4 First Excited State: the Gap 

A better determination of Ac will be for example to compute the energy gap A in the 
spectrum of the Hamiltonian. For values of A near and above Ac the singular behaviour 
of A is characterized by an exponent s, 

A~(A-Ac)^ A>Ac (56) 

which is computed in the disordered region. 

Assuming that the first excited state is translationally invariant and that it belongs 
to the odd sector (R=-l) , we propose the following ansatz, 

where the variational parameter /3 can in principle be different from the ground state 
parameter a appearing in (|52|) . The norm of (|57[) is essentially given by the susceptibility 
at zero magnetic field of a classical Ising model (see appendix D). We have only considered 
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the Id case where the result that we obtain coincides quite surprisingly with the exact 
result, 



Ei{vai, u = l,d = 1) — Eo{vai, u = l,d = 1) 
= Ei{exact, d = 1) - Eo{exact, d = 1) = 2{X - 1) (58) 

For d > 1 the computation of A could be performed using high-temperature expansions 
for the susceptibility. 

4 Concluding Remarks 

We have shown in the study of the ITF that a discrete versions of generalized Jastrow 
wave functions are very adequate to developp a perturbative/variational approach to local 
lattice Hamiltonians. The information encoded in the lowest orders of perturbation theory 
is amplified to all orders by means of a variational method. In this way the ansatzs for 
the ground state of the hamiltonians become less aprioristic, since they are subjected to 
satisfy a set of strong constraints which eliminates much of the arbitrariness. 

An expected weakness of our method shows up for Hamiltonians near criticality, as 
for the energy and magnetization do not exhibit the desired singular behaviour. This can 
be understood from the fact that we do not take into account all the scales involved near 
the critical point. To make further progress one should fertilize the P/V method with 
renormalization group ideas, in the spirit of references |2^, ^ . 
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Appendix ArVariational Calculations in the Ordered 
Phase 



The basic matrix elements needed to compute the energy of the trial wave function (|iO|) 
are the following, 

<^(i)|^(i)>= cosh^/i 
< V"^^^! Enf^n 1^^^^ >= i^cosh^/i tanh/i 

<V'^'^|E„C^nlV'^'^>= Lcosh^-l/l ^ ' 

< ¥'^En,,cr^ ^^+,\^^'^ >= Ld (cosh h)^~^ 
We have assumed that h is real. Using these equations the density energy of ^/''^^■* 
becomes, 

e(var, p = l,d) = — [X tanh h H ni" I (60) 

\ cosh h ) 

The minimum of f|60D is obtained for 



(61) 



\ = 2d tanh/i if A < 2d 
1 = tanh/i if A > 2d 

Introducing (^) into (|60D we obtain (|41|). 

The relevant matrix elements of the trial wave function •\\)'^^^ are, 

<^(')|E„,,^f <J^('^ >= Lrfcosha4-f (a,/i) 



(62) 



i,/^ n n+/i I 

where 2'i(a, /i) is the partition function of a classical statistical Ising model defined in a 
d-dimensional hypercubic lattice. 



{su-,sl} \ n J 



(63) 



{ 

where = ±- 

Zi!^i\a, h) is the same as ZL{a,h) except for the removal of one lattice variable 
located at a single site, say L, together with all possible couplings among this variable 
and the ones in its neighbour. 
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1 



exp a ^ SnSn+/i + hY^ 



(64) 



Finally z{|'°2'^''(a, /;,) is the same as Zi^a^h) except for the removal of two lattice 
variables, say L and L-1 forming a single bond, together with all possible couplings 
among this bond and the remaining lattice variables, namely 



y{bond) 
'L-2 ' 



a, h) 



E 



exp a 



E 



{si,--,SL_2} 



n,n+/i7^L,L— 1 



+ h J2 -^r 



(65) 



To derive eqs.(0) we use the following general formula which relates matrix elements 
of Pauli matrices and statistical sums of Ising models. 



<Tl/(^f,^: 



2 5 



1 



.^A/JlT>=^ E /(Si,S2,---,Sm) 

Sl,---,SM 



where f{xi, ■ ■ ■ , xm) is a generic function of M variables. 
Using eqs(|62D the energy of ip^'^'^ becomes , 



d 



y{bond) 



a, h) 



E{vaT, 1/ = 2, rf) = -A -- In ^^(a, h) - d L 

oh ZL{a,h) 



and the magnetization reads. 



(66) 



(67) 



<a^ > = <^(2)|^Z|^(2)> 



Ziia, h) 

In Id these partition functions can be computed exactly in the limit L 



(68) 



CXD, 



where 



Zr(a,/.) 



Z^l^f'\a,h) 



y{ld,bond) 
'L~2 ' 



a, h) 



cosh h + ^+"'" 
cosh h + 



(69) 



1 

2; = - 
2 



e^cosh h + e~°'\ll + e"^" sinh^/i 
Introducing these equations into (0) we obtain the energy in Id, 



(70) 
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e(var, z/ = 2, ci = 1) = -2 + jt ' W + (71) 

s = v^l + ° sinh^/i (72) 

t = s + Vs2 + e4"-l (73) 

Minimization of ([7I|) with respect to t and s gives, 



A 



9t 

t 



3t2 + 2 + V3t2 + 4 
1 



The magnetization in Id is given, using the variables s and t, by 



(74) 
(75) 



< cr 



t - - + -Vl + - 2s t 



(76) 



Plots of the energy ([7I| ) and the magnetization (|76|), using (|7^ and (|75D, are given 
in figures 1 and 2 respectively. Some numerical results are shown in table 1, where we 



have included for comparison the results obtained in reference I^Qj using a variational 
renormalization group method. In this particular case the variational energy obtained 
through the P/V method is better than the one obtained using the RG method. As for 
the magnetization, in the RG method there is a finite value of A at which it vanishes, 
whereas in the P/V method A has to go to oo. 



A 


energy (var) 


energy (RG) 


< >(var) 


< >(RG) 


0.2 


- 1.010 024 99 [525] 


-1.010 000 


0.994 912 [0] 


0.994 93 


0.4 


- 1.040 399 [417] 


-1.040 306 


0.978 56 [44] 


0.978 83 


0.9 


- 1.212 4 [60] 


-1.210 8 


0.859 3 [125] 


0.860 4 


1.0 


- 1.264 9 [732] 


-1.262 9 


0.8 [0] 


non zero 



Table 1: Energy density and magnetization of the Id ITF obtained by the u = 2 
variational method and the RG methods of reference [^. Numbers in brackets are last 
digits of the exact results. 
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Appendix B: Variational Calculations in the Disor- 
dered Phase 



The matrix elements for the ansatz ('E2D are 



< 



where Zl^q) = Zi^a, h 
The energy reads, 



0) and («) = Zn\a, h = 0). 



{site) I 



(77) 



£'(var, u = l,d) 



d 

-— InZi(a) - XL 
oa 



y{site) I 



a) 



In Id one has, 



(7J 



^l'^^ (a) = cosh''" a 
4^_^f^)(a)=cosh^-2a 



(79) 
(80) 



which yields the energy density. 



efvar, v = 1, = T 



tanh a + 



A 



cosh a 



The minimum of this expresion is obtained for. 



(81) 



i = 2tanha if A > 1/2 
1 = tanh a if A < 1/2 



(82) 



Substituting (|82D into (^) one obtains eq. (|53|) . 

To study the 2d disordered case we need the well known results obtained by Onsager 
in his study on the 2d Ising model [^. Calling f2d = —^i^L^ooj^ogZL^a) the free 
energy per site, U2d = df2d/da the internal energy and C2d = —a'^du2d/da the specific 
heat we have. 



f2d = - In cosh(2a) - ^ /c 



tt/2 



1— sin^ ( 
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M2d = -coth(2a){l + |g' K{q)} 



(83) 



C2d = ^coth2(2a) {K{q)-E{q) 

where 



2 sinh 2a 
cosh^(2a) 
,2/ 



g' = 2 tanh^(2a) - 1 
and -ft'(g) and E{q) are the elhptic integrals, 



(84) 



K{q) 



7r/2 



1 — sin^ I 



Tr/2 



1 — q^ sin 



^5) 



The computation of z^^'^'^^^^^ is a bit more comphcated but it can be related to the 
original partition function through two point correlation functions as follows, 



y{2d,site) 



I [2 + - x(5 + x^) < SqS^^ > +x2(< soS2f,, > +2 < SoSf^^+f.^ >)] (86) 



where x = tanha and Sq, s^^, S2^i, s^^+^j classical Ising variables located at the points 
0, /ii, 2/ii, /ii + ^2 respectively. The correlators entering eq. (|86|) has been computed in 
reference |30|| , 



< Sqs^j >= |cotanh(2a;) 1 + ^q'K{q) 
< SoS^,+^, >= icotanh2(2a) [E{q) + q'K{q)] (87) 
< ^0^2^. >= icotanh2(2a) - {^^' [E\q) - 2q'K{q)E{q) + {q'fK^q)] 

Introducing equations (p3D,(p6|) and ( p?] ) into (^) and minimizing this energy numer- 
ically we obtain figure 3, where we have also include the mean field result (|4l|) for 2d and 
the lower bound given by eq.(p9|). 
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Appendix C: Other P/V solutions 



We mentioned in section 2 the non uniqueness of the solution of the perturbative equations 
(p!0|). This in turn lead us to look for solutions which in the limit A ^ oo will flow to 



the ground state of Hi. This is the reason to impose conditions (pSD- To get further 
motivation for the need of these conditions it is illustrative to investigate another solutions 
to equations (p!0|) than those given in section 3. This will be done only at the lowest order 
u=l. 

In the ordered phase an antihermitean solution of eq. ([To|) is given by, 



which leads to the following ansatz, 



(89) 



^(1) =exp(-f E„^r)l T> 
= On (cosf I T> +sin| I 

where ^ is a real parameter. The variational energy of this state can be computed as in 
appendix A obtaining, 

e(var, zy = 1, d) = - (A sin6' + d cos^9^ (90) 
Whose mimimun is achieved at 6 given by, 

( ^ \<2d 

sin^o = < 2^ ~ 91 
\ 1 X>2d ^ ^ 

Substituting this result in ( pUD we get the same energy as the one given by eq.(^T]), 
which was obtained with the hermitean solution (|40D. Indeed if we compare eqs. (|40|) and 
(^) we see that the trigonometric variable 6 can be related to the hyperbolic variable h 
by means of the Gudermannian function, 

e = gd{h) = 2arctanh(e'^) - | (92) 



The mean field value computation of reference |]T5| is performed in terms of the variable 



6 whose physical meaning is the semiclassical rotation suffered by the spin under the action 
of the external field A, 
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< >= COS0 (93) 

Let us move on to the disordered case. We shall make the following choice of the 
operator f/i, 



= e2 f/i (94) 

where r = (0,?/, z) is a vector with real components and = 1. There are 4 possible 
solutions of the first order perturbative equation (p!0| ) which we display in table 2 together 
with the corresponding variational energies in the Id case. 



Solutions 


e 


r 


e{var, z/ = 1, = 1) 


I 


1 


(0,0, i) 




( -(A + ^) A>i 
-1 A < i 


II 


-1 


(0,i,0) 


2 

27 


A (9 - A2) + (A^ + 3)3/2 


III 


i 


273' 273) 


1 
2 


A + V1 + A2 




IV 


— i 


^273' 373) 


1 
2 


A + V1 + A2 





Table 2: Variational states and energies for the v = l,d = 1 disordered case 



We observe that only solution I, which is the one used in section 3, goes in the limit 
A ^ to the correct value of the energy of the ordered region. Indeed this solution 
satisfies conditions ( [23| ) while the others do not. Table 2 shows that for the ansatz (p^ ) 
there exists no equivalence between hermitean and antihermitean solutions as it happens 
in the mean field case. 



Appendix D: Calculation of the Mass Gap 

The relevant matrix elements for this state are, 

<^ElE„a^|^E>= L(zi-f(/3)-4-f(/3)) (95) 
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where 



ZLi/3) = ^ZLi/3,h)\,=o (96) 



and similarly for Z^]^^^'^\f3) . The expression of the energy of the state (|57|) is, 

Ei(var, = 1, d) = - In Z,{(3) + {ztf - ztf)^ (97) 

A long but straightforward calculation using eqs(^) yields, 

Ei(var, u = l,d= 1) = -L ( tanh /3 + — ^ ) + 2(A - 1) (98) 

\ cosh p J 

whose minimization gives (^). The value of j3 coincides with that of a given by eg. (|52D . 
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Figure Captions 



Fig.l.- Energy versus A in Id: exact (continuous line), mean field (+++) and v = 2 
variational in the ordered regime (■■■). 

Fig. 2.- Magnetization < > versus A in Id: exact (continuous line), mean field (+++) 
and V = 2 variational in the ordered regime (■■■)• 

Fig. 3.- Energy versus A in 2d: lower bound ( |29D (continuous line), u = 1 ordered phase 
(+++) and z/ = 1 disordered phase (■■■)• 
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Fig. 1 
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Fig. 2 
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Fig. 3 
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